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Abstract. We obtain explicit solutions for the density i/py of the first-time 
T that a one-dimensional Brownian process B reaches the twice, continuously 
differentiable moving boundary / and such that fit) > for all t G R"*". 

We do so by finding the expected value of some functionals of a 3-dimensional 
Bessel bridge X and exploiting its relationship with first-passage time problems 
as pointed out by Kardaras (2007). It turns out that this problem is related to 
Schrodinger's equation with time-dependent linear potential, see Feng (2001). 



1. Introduction 

The main aim of this paper is two- folded, given that / G C^[0, oo), and such 



that f"{t) > 0, for all t > 0: 
(1) we compute 

(1) E 



exp|-^'/"(w)^ndw| 



where X is a a 3-dimensional Bessel bridge starting at level a > [where 
a = /(O)] and reaching at some fixed time s, and 
(2) as a corollary, we derive an explicit solution for the density ipf oi the first 
time T that a one-dimensional Brownian process B reaches the moving 
boundary /. 

The derviation of ipf is an old and well studied problem in probability and 
stochastic processes. Yet, explicit solutions existed only in a limited number of 
special cases. 

The appeal of this problem stems not only from its applications (in finance, 
physics, biology, etc.) but also from its mathematical nature. As it turns out the 
problem may be studied through several different approaches: 

(1) Volterra integral equations: Several authors have derived Volterra integral 
equations in the study of the first-passage time problems [the reader may 
consult for instance Peskir (2002) for a detailed technical and historical 
account on the subject]. Yet, in the case in which the equation is not a 
convolution, or its kernel is not separable, numerical procedures are typ- 
ically employed. Thus the results presented in this work provide explicit 
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solutions to previously unknown integral equations. 

(2) Partial differential equations: Groeneboom (1987), Salminen (1988), and 
Martin-Lof (1998) tackled the problem of the quadratic boundary by solving 
a parabolic differential equation, derived from Girsanov's and Feynman-Kac 
theorems. In this paper not only do we generalize this result (an in fact pro- 
vide an alternative proof for this special case), but also show how boundary 
problems are related to Schrodinger's equation with time-dependent linear 
potential. 

(3) Bessel Bridge: Kardaras (2007) notes that the problem of finding is 
equivalent to the study of Several authors, for instance Revuz and Yor 
(2005), have analyzed ([T]), yet our results are new since we do not treat the 
squared process X'^. 

(4) Montecarlo simulation. 

Roughly speaking, the solution to ([1]), is based upon transforming a Cauchy 
problem (derived from Feynman-Kac's theorem), equations © and (fTU]) into that 
of solving the Schrodinger equation (jl3p with time-dependent linear potential [see 
Feng (2001)]. And then proceed backwards to construct the corresponding Green's 
function. 

In Section[2l after stating the problem we: (a) point out the existing relationship 
between a three-dimensional Bessel bridge process and first-passage time problems. 
(h) We present our main results. Theorem 12.31 and CoroUarv 12. 4i which we will 
prove in the remainder of the paper. 

In Section [3l we state the problem in terms of the Feynman-Kac equation. We 
proceed in solving it by finding particular solutions with respect to the backward 
(t, a) and forward (r, &) variables, in Sections H] and [5] respectively. We conclude in 
Section [6] by using these results to construct the necessary Green's function G. 

Both the backward and forward equations in Sections |4] and [5] are dealt with in 
a similar fashion: (1) Through algebraic transformations we obtain Schrodinger's 
equation with time-dependent linear potential. (2) We solve the corresponding 
Schrodinger's equation, by reducing it into the heat equation. 



Problem 2.1. The main motivation of this work is finding the density LpT of the 
first time T , that a one- dimensional, standard Brownian motion B reaches the 
moving boundary f : 



where f is twice, continuously differentiable, and f"{t) > for all t > 0. 

Alternatively, from Girsanov's theorem, the "problem" may be restated as: 
Proposition 2.2. Given that 



2. The problem 



(2) 



T :=M{t>0\Bt^ fit)} 



(3) 




t > 0, a > 
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is the density of the first time that B reaches the fixed level a e (0,cxd), the process 
X is a 3-dimensional Bessel bridge, which has the following dynamics: 



dXt ^dWt + ( ^ ^ I dt, Xo^a, te [0, s). 

\Xt s-tj 

(Where W is a Wiener process.) Then, the distribution of T equals 



p(r <t) = 







[ E 


cxp 1 


'0 





X exp 



f"{u)Xudu \ 
{f'{u)fdu~f'{0)a\Ms)ds 



Proof. From Girsanov's theorem, and the fact that the boundary / is twice contin- 
uously difFerentiable, the foUowing Radon-Nikodym derivative: 



^{t) exp^-/ f'{s)dB,~^- I {f'{s)fds 



exp { -f'{t)Bt + f"{s)B,ds ~-J^ {f'{s))'ds } , 



is indeed a martingale and induces the following relationship: 



p(r < t) 



= E 



dP, , 



E 



dP , , 
^ T 
dP 



T = s 



ipa{s)ds, 



where the second equality follows from the optional sampling theorem, and the 
third from conditioning with respect to T under P, and tpa is defined in ([3]). 

Next, Kardaras (2007) shows that calculating [see for instance Chapter 11 in 
Revuz and Yor (2005), for a detailed overview of 3-dimensional Bessel bridges] : 



(4) 



E 



^(T) 



T = s 
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is equivalent to finding the expected value of a functional of a 3-dimensional Bessel 
bridge X, which at time starts at a and at time s equals 0. Indeed, 



E 



■E 



exp{-f{T)a+ / f"{u)Budu 



if'iu)fdu 



T = s 



E 



1 



exp <j -f'is)a + I f"{u)ia - Xu)du - ^ / i.f'{u))^du 



exp|-/'(s)a + (/'(s)-/'(0))a- / f"iu)X.,du 
{f'iu)fdu 



(5) 



1 

2 Jo 
exp 



f'iu)X^du 



exp 



{r{u)fdu~f'{0)a 



thus the process a — X equals at i = 0, and at s reaches level a for the first time, 
as required. □ 

We conclude from equation ([5]) that our goal is to compute the following expected 
value: 



(6) 



E 



exp 



f"{u)X^du 



Equation ([6]) has the following solution: 
Theorem 2.3. Let the boundary f e €^[0, oo) satisfy f"{t) > for all t. The 



E 



exp 



f"{u)Xudu 



v{t, a) 



where 

(7) v{t,a) 
given that 



G(t,a;s,b)db, 0<t<s, a e R+\{0} 



G(t, a; r, 6) = — ^H{t,a;T,b), 0<t<T<s 



fais - t) 



where 



H{t,a;T,b) = 
1 

VMr - t) 



exp 



if'{u)rdu~nT)b + f'{t)a 



exp ■ 



{b-a-J[ f{u)duy 



exp • 



{b + a~ J[ f'iu)du) 



2^1 



2{T-t) 

Where (pa is defined in (0)- 

Thus the density ipr of the stopping T in equation ^ equals 



2(r-t) 
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Corollary 2.4. Let B be a one- dimensional standard Brownian motion B and let 
the stopping time T be as in Then, the density ipx equals: 

^t{s) = ^-(O, a)^a{s) exp |-i ^ {.nu)fdu - /'(0)a| . 

where v and if a are defined in ^ and ^ respectively. 



3. The Feyman-Kac equation 



As the reader may suspect from equation ([6|) and the Green's function appear- 
ing in equation ([7]), in Theorem 12.31 we will follow a P.D.E approach to solve our 
problem. To do so, we will first introduce Feynman-Kac's theorem and some neces- 
sary tools, and will conclude this section by introducing the corresponding Cauchy 
problem (see Problem 13. 3p . 



b{t,a) — {t,a) 



If we let the operator At satisfy 
{Atv){t, a) 



1 2/ N^^^^/ N 

-a {t,a) — it,a) 



2d^ 



{t, a) + 



1 



dv 



— t J da 



{t,a) 0<t<s, ae M+\{0}, 



Then, we may introduce a simplified version of the celebrated Feynman-Kac theo- 
rem which relates parabolic partial differential equations with stochastic differential 
equations [for a detailed and general overview of this and related topics the reader 
may consult for instance: Chapter 5, in Karatzas and Shreve (1998)]: 



Theorem 3.1 (Feynman-Kac). Given the Gauchy problem 
dv 

- +kv = AtV + g; [0, T) x M 
ot 

v{s,a)^f{a)- aeM+\{0} 



where k{t,x) : [0, s] x R+ — > [0, oo), and functions f and g satisfy the standard 
assumptions [see Theorem 5.7.6 in Karatzas and Shreve (1998)], and v satisfies 
the polynomial growth condition. Then, the function v{t,a) admits the .stochastic 
representation 



f{Xs) exp <^ - / k{u, Xu)du 



k(u, Xy)du > dv 



v(t,a) = E* 



g[v,X^) exp 

Jt 

on [0, s] X M+; in particular, such a solution is unique. 

The following definition, will guide us on the construction of our Green's function 

G: 

Definition 3.2. [Definition 5.7.9 in Karatzas and Shreve (1998)] A fundamental 
solution of the second-order partial differential equation 



(8) 



i5u , , 

- — + few = Atu 
at 
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is a nonnegative function G{t,a;T,b) defined for 0<t<T<s, aCz K+, b € 
with the property that for every f £ Co(M+), r G (0, s], the function 



u{t, a) 



G{t,a;T,b)f{b)db; < t < r, a G M+ 



IS 



bounded, of class C^'^, satisfies and 

\imu{t,a) ^ f{a), a G IR+. 

For fixed (r, b) G (0, s] x R+, the function 

ip{t, a) := G(i, a; r, 6) 

is of class C^'^ and satisfies the backward Kolomogorov equation in the backward 
variables (t, a). And, for fixed (t, a) G [0, s) x M+ the function 

Tjj{T,b) := G(t,a;T,6) 

is of class C^''^ and satisfies the forward Kolmogorov equation (also known as the 
Fokker-Planck equation) . 

Hence, from equation ([5]) it follows that: 



E' 



0,a 



exp 



f"[u)Xudu 



= w(0,a). 



where 



/(O)exp - / k{u,Xu)d 



v{t,a) := E* 



for k{u,a) := f"{u) ■ a, and u(T, a) = /(a) = 1, i.e., 

Problem 3.3 (Bessel's Cauchy problem) . Equation (0) satisfies the following Cauchy 
problem: 



(9) 



(10) 



^{t,a)+r{t)avit,a) = ^£(i,a) 

1 a \ , , 
a s — t J oa 

v{s,a) = 1; 



m [0,s) X R+\{0} 



a G 



A{0}. 



Definition 13.21 hints us towards the path we should follow. We should first find 
a general solution to the backward-Kolmogorov equation, in Section [4| 



1 0"^ 

-^(t,a) + f"{t)a^{t,a) - i^^H^a) 



+ 



1 



a \ dip 



(t, a); 



^ a s — t J da 
and to the forward-Kolmogorov equation, in Section [5j 

(r, b) + 4 (T)bV(T, b) = 2 



[0, s) X R+\{0} 



96 



1 



i/;(r, 6) 



in (0,s] X M+\{0}. 
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Which after proper change's of variable and algebraic transformations will equal 
respectively: 

and y and z are function of a and & respectively (and will be determined later). 
Since both Lp and V' may be reduced to the Heat equation, and the variables a, h 
are only defined on the half line. Then the Green's function G which relates these 
two solutions will be of the form: 

'^■"--'^^ (-{-It^} --{-^}) 

Thus, the construction of the particular solution G will just be a backwards proce- 
dure. 

4. Schrodinger's backward equation 

In this section we obtain a particular solution to the backward equation ([9]) 
by first transforming it into Schrodinger's backward-equation with time-dependent 
linear potential through algebraic transformations. 

Proposition 4.1. Equation (0i satisfies the following relationship 

(11) ifit, a) = if' {t, a) {A{t) exp [B{t)a'] } , 
where 

^^^^ ^ W^) 

A{t) = c • (s - tf/^ 

(12) a/"(t)/(t, a) - \ip\,{t, a) - -^l{t, a) - ^\{t, a) = 

2 a 

Proof . Let us first start by analyzing the transformation in equation (111^ : 

|^(t,a) = ^l{t,a)A{t)c-^Y>[B{t)a^] 

'2aLp\t,a)A{t)B{t) exp [B{t)a'^] 



Q^it,a) = ipUt,a)A{t)e^p[B{t)a^] 

+2aipl{t,a)A{t)B{t) exp [B{t)a'^] 
+2(p\t, a)A{t)B{t) exp [B{t)a'^] 
+2aipl{t,a)A{t)B{t) exp [B{t)a'^] 
-|-4aV^(i, a)A{t)B^{t) exp [B{t)a'^ 

^{t,a) = ip\t,a)At{t)e,cp[B{t)a^ 



dt 



(t,a)A(t)Bt(t) exp [B(i)a^ 
-|-(pt(i,a)A(t) exp [B{t)a^ 
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Substituting these equations into ^ we get 



.__ _ „ B,{t) ^-+af {t)^-— + a-B{t) 

+B{t) + a^B{t) + 2a'B\t) +(-- 4 

- 2aB{t), 
a s — t J 



equating to zero we have 



_ a B,{t) - _ + a/ (i) - --^ - 2a-,Bit) 
-B(t)^2a'B'{t) ■ " 



a Lp^ s — t ip^ 



a" 



~2B(t) + 2 B(t) = 

s — t 



finally, we factorize 



+ 3B{t)) ~ (Bt{t) + 2B\t) - 2^—B{t) 



-a^ { 2B{t) — ) + af"{t) - - -^^ _ ri ^ q. 

tp^ \ ^ ' s-tj ^ ' 2 (^1 atp^ (^1 

This equality holds for every a > if and only if 

B{t) 



2is - t) 
A{t) = c • (s - tfl^ 

af"{t)ip\t, a) - \ipl,{t, a) - \l{t, a) - p>\{t, a) - 
2 a 



as claimed. □ 



4.1. The Schrodinger's backward equation. If we let pP = a ■ ip^, then 



1 1 2 

a 

1 1 2 1 2 

fa = -Pa - —'P 

a 

1 "^2 -^2'^2"^2 

Vaa = -P>aa~ —'Pa + —"P - —"Pa, 

a a-^ 



af"{t)v^{t, a) - ]^vla{i^ a) - Lpl{t, a) = 0. 



after substituting in equation (|17p 
or, equivalently 

(13) - = ^^(^.«) - ^f"{t)v\t,o) 

which is Schrodinger's backward-equation with time-dependent linear potential [see 
for instance, Feng (2001)]. 
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Theorem 4.2. Given the backward Kolmogorov equation in \13\l: 
(14) 



~^-^{t,a) = \^{t,a)-af"(t)^\t,a) 



we have the following fundamental solution, 

{a- fit)? 



(15) ^^(t,a) = 



exp 



2t 



-\ I {nu)fdu + nt)a 



where 3 



-1. 



Proof of JJ5)) . First set ip^{t,a) — X(t,a)e^''*^°', where P will be determined later 
and compute 



da 
5a2 



it,a) 
{t,a) 



aA(OA(t,a)e^W'^ + |^(t,a)e^(*)» 

mKt,ci}e^'^'^^ + ^{t,a)e^^*')- 
oa 

f^'{t)X{t,a)e^^'^- + m^it,a)e^^*^'^ 



Ik? 



{t, a)e' 



I3(t)a 



dX 

da 
dX 



,l3{t)a 



(3^{t)X{s, a)e^W'^ + 2/3(i)^(t, a)e''W" 

oa 



next, substitute into equation (|T4I 
dX 



dt 



\ d^X 1 dX 

{t,a) = - — (i,a) + -/32(t)A(t,a)+/3(0^(t,a) 



+a{Pt{t)-f"{t))X{t,a) 
Next, if we let y = a — v{t) and set A(i, a) =^ u{t, y), then 



to delete the term 



set vt{t) = [3[t), hence 



1 O^U 1 C^ti 

+ iy + vit)){Ptit)~nt))Xit,y) 



du , 



du , , 



1 



+ (y + «(t))(A(i)-/"(0)A(t,y) 



furthermore, setting /3t(t) = /"(i) we get 



du , , 1 c^^u , , 1 ^9/ N / X 
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which is a parabohc differential equation which now may be transformed into the 
heat equation by letting G{u) — 1/2/3^ (u) and 



it follows that 



which has general solution equal to: 



exp 



I.e. 



g " 2 /o f 

9 



^ p-|/o(/'(«))'*V'(i,y) 
: exp 



and 



X{t, a) = 



/27rt 



/27rf 



/27rt 



exp 



exp 



7/2 1 /■* 



2t 2 7o 



(a-t^_ 1 
2t 2 

2t 2 



{f'{u)fdu 



inu)fdu 



Hence 



/27rf 



/27rf 



exp 



exp 



2t 



-2 / (/'H)'d« + /3Wa 



2t 



-t I {f'{u)fdu + nt)a 



5. Schrodinger's Forward equation 
From equation ^ we will now study its dual, the Fokker-Planck equation: 



dr 



1 9V 

2 952" 
_d_ 
~db 



b s ~ T 



^{r,b) 



(The techniques used in this section are equivalent to those used when treating the 
backward-equation, with the obvious modifications.) That is 

(16) (r, h) - T^7Jiy(T,6) 



2 962 
'1 



6 s — T J db 
1 1 

62 S — T 



/^'(r)6W'(r,6) 



ON SCHRODINGER'S EQUATION, 3-DIMENSIONAL BESSEL BRIDGES 11 

Proposition 5.1. Equation il6\) satisfies the following relationship 
(17) b) = V'(t, b) {A{t) exp [B(t)&2] } , 

where 

B{r) ^ 



2{s - t) 
A{t) =c-(s-r)-i/2 

(bf"{T) - 1 - Vi(r, 6) - iv6\(r, b) + i0i(r, 6) 

(18) +V4(^,fe)-0 

Proof . Let us first start by analyzing the transformation in equation (|18p: 

^(t,6) = V6'(r,6)A(r)exp[i?(r)62] 

+26?/'\t, 6)A(t)B(t) exp [B(t)&2] 
^(r,6) = 0i,(r,6)A(r)exp[i?(r)62] 



+26^6 (t, 5)A(t)B(t) exp (^(t)^^) 
+2V'\r, 6)A(t)B(t) exp [^(t)^^] 
+26V'b A(t)S(t) exp [5(t)6^] 
+4&V^(t, b)A{T)B^{T) exp [B(r)&2] 

|^(t,&) = V^i(r,6)A.(r)exp[S(r)&2] 

+6V^(t, 6)A(t)B^(t) exp [B(t)62] 
+V4(T,fe)A(r)exp[B(T)62] . 



Substituting these equations into (|16p we get 
^r(r) 



A(r) 



+ 6|ii?(r) + B{r) + b^B{r) + 2b- B\r) 



1 5 \ V'b /I 



6 s — T J ^j^ \b 



2bB{T), 



equating to zero we have 



A{t) ' V^i V 



1 V'hfc 

2 










S — T J 



1 6 



6 s — T J 



] 2hB{T) = 0, 
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finally, we factorize 



^^^^^ ' Bir)) +b^( Br{T) - 2B\t) - 2-^B{t) 



A{t) ' 7 V -5-^ 



This equality holds for every 6 > if and only if 

B{t) = 



2{s-r) 
A{t) =c-(s-t)-i/2 

(bf"{r) - 1 - ^) V^i(r, 6) - \i^l{T, b) + \i^l{T, b) 

+V;^(r,a) = 



as claimed. □ 



5.1. The Schrodinger's forward equation. We first start with the following 
transformation -0^ = b ■ ip'^, i.e. 







^bb 


= 2^1 + bi^l^ 







which after subtitution in ()18p equals 

s — r 2 

S — T I 

(19) f6/"(r)-^')^'"iV',\ + V'? - 0. 



2 



Next, setting t/j^ equal to: 



we have that 



s — r 
"066 = -^-^Ib- 

S — T 
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Thus subtituting into (fTS]) 



bf'ir) J2 1 72 1 1 72 ^ 1 72 ^ 1 72 n 

s — T (s — r)^ 2 s — r (s — r)^ s — r 

bf'ir) 72 1 1 72 1 72 

S — T 2 S — T S — T 

(20) ^/"W'-^^b\ + V^? = 
which is Schrodinger's forward equation. 
Theorem 5.2. Given the forward equation as in Ii20\) 

flip 1 r)2,/'.2 

(21) ^{r,b) = i^(r,6) - 6/"(r)^2(r,6) 

we have the following fundamental solution 

(22) 4''{r,b) = -^exp|-i^::|^^ + i^^(/'(u))2du-/'(r)6 

Proof of IM)) . First set ip'^{T,b) — \{T,b)e^^'''^^'' , where /3 will be determined later 
and compute 

^(r,6) = -6/3.(r)A(r,6)e-''M^ + |^(r,6)e-'^M^ 
^(r,6) = _;3(r)A(r,6)e-^M'' + |^(r,6)e-^W'' 



^(r,6) = /3^(r)A(r,6)e-^(^)^-/3(r)|^( 



+ 0(r,5)e-^(^)^-/3(r)|e-^(^)^ 
/32(r)A(r,6)e-'^M^-2/3(r)|^(r,5)e-^M'' 
+ |^(r,5)e-'^(^)^ 



next, substitute into equation (|2ip 

+6(/3r(T)-/"(r))A(T,6) 
Next, if we let z = — i'(r) and set A(r, 6) — u{r^ z), then 
(9w du 1 d^u 1 

— (t,z)-v,(t) — (t,z) = __(t,z) + -/32(t)m(t,z)-/3(t) — (t,z) 

+ (z + «(r))(/3,(T)-/"(T))zi(T,z) 

to delete the term 
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set Wt(t) — (}{t), hence 

+ iz + v{Tmr{r)-f"{T))u{T,z) 

furthermore, setting Pt{t) — /"(t) we get 
du , 



dr 



which is a parabohc differential equation which now may be transformed into the 
heat equation by letting G{u) = 1/2/3^ (u) and 



it follows that 



whose solution is 



/27rr 



cxpi-- 



I.e., 



and 



A(r,5) 



1 



/27rr 
1 

/27rr 



/27rT 

exp 
exp 



2r 2 

ib-f{r)r 1 r... 

2r 



Hence 



1 



/27rT 
1 

/27rr 



exp 
exp 



+2/ U W) du - l3{T)h 



6. General solution 

From the backward (fT5|) and forward (|22|) solutions of Schrodinger's equation, 
and the fact that a, 6 £ M+\{0} we have 

Proposition 6.1. Schrodinger's equation, see equations {1^^ and k21\). with time- 
dependent linear potential and such that a,b G M"'"\{0} has the following Green's 
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function 

(23) H{t,a;T,b) 



^ expji / if'{u)fdu-riT)b + f'{t)a 



{ {h-a-ll r{u)duy \ f {b + a-J[ f'iu)du) 

exp < - ^ — — ^ 7 — > — exp 



2 "\ H 



2(r-t) j ^1 2(r-t) 

where < t < t < s. 

Proof. Recall that the backward and forward equations may be reduced corre- 
spondingly to: 

where y ^ a - f{t), and z = b- f{T), for < t < r < s, and a,b £ R+\{0}. Thus, 
the Green's function G that links these two solutions should be of the form: 

Finally, equation (|23p follows from equations ([M]) and pT|) . 
Which alternatively leads to 

Proposition 6.2. Oitr initial equation (0) /las f/ie following Green's function: 
G{t,a;TM^ '^'^i^^l H{t,a-T,b) 

•fa^S - t) 

Proof . It follows from equations pTjl . (fT7|) . and and verification of Definition 

Proof of Theorem 12.31 It follows from Proposition [621 
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